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Abstract. We consider gravity theory with varying speed of light and varying gravitational 
constant. Both constants are represented by non-minimally coupled scalar fields. We ex¬ 
amine the cosmological evolution in the near curvature singularity regime. We find that at 
the curvature singularity the speed of light goes to infinity while the gravitational constant 
vanishes. This corresponds to the Newton’s Mechanics limit represented by one of the vertex 
of the Bronshtein-Zehuanov-Okun cube [1, 2]. The cosmological evolution includes both the 
pre-big-bang and post-big-bang phases separated by the curvature singularity. We also inves¬ 
tigate the quantum counterpart of the considered theory and find the probability of transition 
of the universe from the collapsing pre-big-bang phase to the expanding post-big-bang phase. 
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1 Introduction 

Among the most popular theories admitting varying physical constants are those involving 
the variable fine structure constant a, the speed of light c, and the gravitational constant G. 
A lot of important motivation for such theories have been provided on both theoretical and 
experimental basis [3]. 

The most important framework for the variability of the gravitational constant is pro¬ 
vided by the scalar-tensor theories, where the effective gravitational constant is determined 
by the local value of the scalar field [4], 

Theories with variable speed of light can generally be divided into two groups. The first 
group contains theories that assume the dependence of the speed of light on the space-time 
coordinates. Such an assumption unavoidably leads to violation of the Lorentz invariance [5- 
10]. In the second group there are theories where the speed of light depends on the energy 
scale [11]. The most radical example of a theory violating Lorentz invariance assumes the 
existence of a preferred frame usually identified with the cosmological frame [5, 6]. In these 
theories the dynamics of the gravitational field, matter fields and the scalar field representing 
the speed of light is described by gravity theory with minimal coupling in the preferred frame. 
It is also assumed that the scalar field which determines the value of the speed of light does 
not contribute to the energy momentum tensor. Consequently, in the preferred frame the 
gravitational field is described by the standard Einstein equations with the speed of light 
being merely replaced by the above mentioned scalar field. The bimetric VSL theory [7-9] 
is another example of a theory that brakes Lorentz invariance. In this framework the two 
separate metrics are introduced. One of them describes the geometry of space-time while the 
other couples to the matter. As a consequence, in such a theory the speed of gravitational 
waves may differ from the speed of massless particles. Another interesting subclass of theories 
with varying speed of light exploits a concept of local Lorentz invariance which is applicable 
when the speed of light is allowed to vary [10]. The second class contains theories assuming 
a modified dispersion relation, wherein the modification becomes important for the energy 
scales comparable with the Planck scale. As a consequence in the considered model the group 
velocity of light becomes dependent on the energy scale. Originally, such theories assumed 
the modified dispersion relation to be fulfilled in the preferred frame and hence they explicitly 
violated the Lorentz invariance [12-14]. As it was shown in [15, 16], it is possible to preserve 
the Lorentz invariance (its non-linear form), while still assuming modified dispersion relation. 

In refs. [17-19] the tunneling from nothing to a de Sitter phase within the quantum 
cosmology with varying physical c and (or) G constants was discussed. In ref. [20] the theory 
of gravity with varying G and c, with the minimal coupling in the cosmological frame was 
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considered. Neither c nor G did contributed to the energy momentum tensor. Consequently, 
the cosmological equations were of the form of the standard Friedmann equations with both 
constants replaced by arbitrary time dependent functions. By using the quantized theory 
described by the one-dimensional Wheeler-De Witt equation the probability of tunneling of 
the universe from nothing was derived. A similar approach was adopted in ref. [21] where it 
was shown that the cosmological constant can be considered as an eigenvalue of an appropriate 
Sturm-Liouville problem of a varying speed of light theory. We think it would be advisable to 
consider a theory where both varying c and G represent additional degrees of freedom. This 
could be achieved by replacing c and G with two separate scalar fields. 

Our paper is organized as follows. In Sec. II we formulate a theory with varying 
c and G where both constants are represented by non-minimally coupled scalar fields and 
find the classical solution near the curvature singularity. We show that the cosmological 
evolution is equivalent to the process of the one-dimensional scattering of the point particle 
on the exponential potential barrier. In Sec. Ill we solve the Wheeler-De Witt equation 
for the quantum counterpart of the considered theory and find that the quantum behavior 
is equivalent to the one-dimensional quantum stationary scattering of the point particle on 
the exponential potential barrier. By calculating the probability of scattering we find the 
probability of the transition from the pre-big-bang collapsing phase to the post-big-bang 
expanding phase. In Sec. IV we give our conclusions. 


2 Non-minimally coupled classical varying c and G theory 


We consider a theory which admits varying speed of light c and varying gravitational constant 
G. This is achieved by connecting both constants c and G with the two scalar fields (p and 
■ip. Additionally, we include separate kinetic terms for each scalar field. Such a theory can be 
represented by the following action: 


S = J^/^(^^'^[R + A + uJ{^|,(p^^cp + ^|,'^p^^"'^p)]d^x, ( 2 . 1 ) 

where a; is a constant parameter and the relations between varying c(x^) and G{x^) and the 
scalar fields <p{x^) and p){x^) are given by = e'^ and G = . By the application of the 

following field transformations: 
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the action (2.1) can be rewritten in the form of the Brans-Dicke action: 


( 2 . 2 ) 




5{R + A) + + 53^158^^15 


d^x , 


(2.3) 


where the constant parameter uj can now be identified with the Brans-Dicke parameter. 
We consider a Friedmann universe described by the FLRW metric: 


ds^ = —N{tY{dx^’Y + a{t) 


dr"^ 


+ Ydn^ 
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where N{t) is the lapse function, a{t) the scale factor and fc = 0, ±1. By inserting the metric 
(2.4) into the action (2.3) and then integrating by parts we get: 


S=^ / ( ~a'S' - ^aa'^ + ka6N + A6a^N 


Sir 
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(2.5) 


where ( )' = By choosing the lapse function as N = a^5 and further by application of 
the helds transformations given by: 


X = ln(a\/5), 


y = — Inh. 

2A 


( 2 . 6 ) 


where A = 


^ „ , one arrives in the following form of the action: 

Y 1— lu) ‘ 


Svr ./ 


+ y'2 ^ ^ ^^2{X-AY)\ _ 


d/2 


(2.7) 


From now on, we restrict our considerations to the spatially flat Friedmann universes and put 
k = 0. Further helds redehnition I = AY — 3X, J = 3y — AX and B = V^/3 gives: 


5 = j dx^ [m( J'2 - /'2) + Ae-2^ - , 

where we have dehned the following constant parameters: 

3^0 Vh 


( 2 . 8 ) 


V^) = 


m = 


A = VoA. 
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Equivalently, our theory can be described by the following Hamiltonian: 




1 


- -I) - -I 


m 


— Ae 


-21 


(2.9) 


where irj = 2mJ', irj = —2m/', and itb = —2B' are the canonically conjugated momenta. 
The form of the hamiltonian (2.9) clearly shows that irj and itb are conserved. Consequently, 
the classical evolution near the singularity can be reduced to the problem of motion of the 
point particle in the one-dimensional exponential potential. The associated Hamilton equa¬ 
tions read as: 


B' = 


/' = 


dH 

dlTB 

dH 

diri 


'2m ’ 


H - n 

/ dH - .j 

., = -^ = -2Ae- 


( 2 . 10 ) 


jt ^ ^ ^ ^ 

dir j 2m ’ 


= 


dH 


= 0 . 


The solution of the Hamilton equations (2.10) is given by: 

B = -^x^ + P, 

/ = In sinh | y/[A^ — 9)Aj;‘^ | , 


( 2 . 11 ) 













where E and P are new integration constants. This solution expressed in terms of the scale 
factor a and the held 5 is: 


1 

sinh'^ 1^ (^_^2 _ 9)AxO 
sinh^ I {A^ — 9)Aa:° | 

where 


( 2 . 12 ) 


D = 
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and > 9. It is also possible to relate variable with the proper time r encountered by 
the comoving observers since both quantities are related by: 


dr 


V—ds'^ 

c(x°) 


(2.14) 


where ds'^ = —N‘^{dx^)‘^. The relation (2.14) is a postulated generalisation of the general 
relativistic dehnition of the proper time for the case of varying speed of light c and it is assumed 
to be fulhlled in the cosmological frame. A justihcation for this particular generalisation of the 
proper time dehnition arises from the fact that we require light rays to follow null geodesics. 
The different way of introducing varying c into the metric was presented in ref. [22] where it 
was argued that varying c merely rehects the rescaling in the time coordinate and thus does 
not affect the value of the physical speed of light. In the dehnition (2.14), however, c does 
represent the physical value of the speed of light. Noticing that: 


sinh I {A^ — 9)Kx^\ 

we obtain: 

arctanh (, for < 0 , 

v/(A2 - 9)A V / 

x^ = — arctanh (, for > 0 . 

V^(A2 - 9)A V / 


(2.15) 


(2.16) 


where we have dehned dx^ = c{x^)dT. In Fig.(l) we have plotted the behaviour of the 
scale factor a (gray), the speed of light c (black), and the gravitational constant G (dashed) 
for the pre-big-bang x® < 0 and post-big-bang > 0 phases. At x^ = 0 an interesting 
behaviour occurs: as the universe approaches the singularity (a —?■ 0), the gravitational 
constant G goes to zero while the speed of light c reaches inhnity. This corresponds to 
the Newton’s Mechanics limit represented by one of the vertex of the Bronshtein-Zelmanov- 
Okun cube [1, 2]. An analogous behavior is encountered in the case of the ekpyrotic [23, 24] 
and cyclic scenarios [25, 26], where in the high curvature regime the universe, after having 
undergone the phase of the accelerated contraction, faces the big-crunch/bang singularity 
and then bounces to enter the standard expanding phase. On the other hand, we encounter 
a qualitatively different behavior in the pre-big bang scenario within the framework of string 
cosmology where the pre-big bang phase corresponds to expanding and accelerating universe 
approaching the strong coupling regime [27]. 
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Figure 1. Qualitative behaviour of a (gray), c (black) and G (dashed) before (a;° < 0) and after 
(a;° > 0) the big-bang singularity. 


3 Quantum mechanical behavior near the curvature singularity 


By examining the classical solution we find that the high-curvature (near big-bang) regime 
appears for I —)• oo, which is equivalent to —>• oo for —)• 0, while the low-curvature 
(far away from big-bang) regime appears for I —)• —oo, which is equivalent to —)• 0. 

The latter regime corresponds to —>■ — oo or to —)> oo for the pre and post-big-bang 

branch, respectively. Both regimes are characterized by the following asymptotic values of 
the momentum ttj: 


TTl 




collapsing pre-big-bang 
expanding post-big-bang 


in the high-curvature limit and 


TTl 




collapsing pre-big-bang 
expanding post-big-bang 


in the low-curvature limit. In order to quantize the considered theory we make the following 
substitution: 


d ^ d 

TTJ vrj = , TT/ -)■ TT/ = -z — 

oJ ol 




TTb = —I 
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(3.1) 


The resulting Wheeler-De Witt equation is: 
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(3.2) 


It is interesting to point out that the sequence of the helds transformations given by (2.2) 
and (2.6) reduces the problem of the factor ordering. We look for the solution of eq. (3.2) 
which is separable [28]: 

d> = a(J)7(B)/3(/), 
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where 


and 


i2k 


Vo 7 


a{J) = e 

7(B) = 

I3{I) = J-i^^ {7rfe~^) 


(3.4) 


FoA = k^ + l\ 


(3.5) 


The wave function 4>, asymptotically for I —>■ oo, is an eigenfunction of the momentum 
operator tt/ since: 

717/3(1) = 777-/3(1) . (3.6) 

Thus, in the high-curvature regime <1> contains only those modes which represent collapsing 
universe [29, 30]. In the low curvature limit (/ —)• —oo) the Bessel function Jin can be 
decomposed in the following way: 


where 


and 


Jin{z) = 4'i + Ta , 



TT/^'l = —zTl 
TT/'I'a = z'^2 


(3.7) 


(3.8) 


(3.9a) 

(3.9b) 


where n = —7r|° and z = 7r|°e ^ a^- 9 ^ ^ (3.9a) and (3.9b) therefore imply 

that Ta and Ti correspond, to the pre- and post-big-bang branches in their low curvature 
evolution phases, respectively. The probability of the transition form the pre-big-bang low 
curvature initial state to the post-big-bang low curvature hnal state is given by the reflection 
coefficient: 

R = . (3.10) 

On the contrary, in the scenario provided by the string quantum cosmology [30, 31], the 
initial state from which the quantum transition to the expanding phase occurs, represents an 
accelerating expanding universe. 
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4 Conclusions 


We have considered a theory with varying c and G where both constants are represented by 
non-minimally coupled scalar fields. Such theory is equivalent to the Brans-Dicke theory. We 
have examined the cosmological evolution near the curvature singularity. It turned out that 
dynamics of the expansion can be reduced to the problem of one dimensional scattering of the 
point particle on the potential described with an exponential function. We have found that 
the cosmological evolution contains both the collapsing pre-big-bang phase and the expanding 
post-big-bang phase with the both phases separated by the curvature singularity. At the 
curvature singularity c reaches infinity while G teds to zero. Thus, we have shown that the near 
curvature singularity evolution phase takes place well within the Newton’s Mechanics regime 
given by G = 0 and c = oo. We have also examined the quantum mechanical behavior near 
the curvature singularity by examining the Wheeler-De Witt equation. We have shown this 
to be equivalent to the quantum stationary scattering of the point particle on the exponential 
potential barrier. The solution of the Wheeler-De Witt in the low curvature regime is the 
superposition of the wave function that represents a particle with positive momentum moving 
towards the potential barrier and the wave function that represents a particle scattered on the 
barrier. The former correspond to a collapsing universe in the pre-big-bang phase while the 
latter represents an expanding universe in the post-big-bang phase. We have also calculated 
the probability of scattering, which is equivalent to the probability of the transition from the 
pre-big-bang phase to post-big-bang phase. 
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